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!•  Introduction 

The  linear  theory  of  elasticity  of  anisotropic  media  has  long  been 
the  sxibject  of  ntimerous  investigations.  In  particulari|  the  problems  of 
torsion  and  flexure  of  beams  for  various  underlying  anisotropic  stress- 
strain  relations)  as  well  as  the  plane  problem  for  an  orthotropic 
mediian,  have  been  treated  extensively  in  the  literature  which  now  contains 
a wealth  of  significant  solutions  of  such  problems.^ 

The  abandonment  of  the  assumption  of  elastic  isotropy  in  three- 
dimensional  problems)  on  the  other  hand)  leads  to  greater  analytical  com- 
plications. Here)  apparently)  only  the  case  ii>  which  the  material 
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possesses  an  axis  of  elastic  syianetry  ("transverse  isotropy"  )»  as 
illustrated  by  crystals  of  the  hexagonal  system)  has  received  successfiil 
attention.  The  first  izivestigatlon  of  this  kind  is  due  to  Hlchell  H 
(1900))  who  obtained  the  solution  corresponding  to  a transvez^ely  iso- 
tropic semi-infinite  elastic  body  under  arbitrarily  prescribed  surface 
tractions  on  its  plane  boundary. 

The  results  comsiunicated  in  this  paper  were  obtained  in  the  course 
of  an  investigation  conducted  \mder  Contract  M7onr-32906  with  tJie  Office 
of  Naval  Research)  Department  of  the  Navy)  Washington  2$)  D. 

\ bibliograjiiy  of  two-dimensional  investigations  based  on  aniso- 
tropic stress-strain  lawS)  is  beyond  the  scope  of  the  present  paper. 

^his  terminology  is  used  by  Love  Cil  ) p.  l60.  Numbers  in  brackets 
refer  to  the  bibliography  at  the  end  of  the  paper. 
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Further  progress  in  this  connection  was  made  by  Lekhnitsky  ^3^  (I9I4O), 

irtjo  stated  a stress-function  approach  to  problems  characterized  by 
torsionless  axisymaetry  for  the  case  in  which  the  axis  of  elastic  synnetry 
coincides  with  the  axis  of  stress  symmetry.  Lekhnitsky 's  solution  of  the 
displacement  equations  of  equilibrium  in  terms  of  a single  stress  function 
satisfying  a fourth-order  partial  differential  equation,  is  a generaliza- 
tion of  Love's  approach^  to  isotropic  rotationally  symmetric  pf*oblems, 
which  in  turn  may  be  regarded^  as  a specialization  of  Galerkln's  M 
general  solution  of  the  three-dimensicnal  field  equations.  Lekhnitsky 
[[3]  employed  his  approach  in  sectoring  several  technically  Important  paor- 
ticular  solutions. 

Subsequently,  Ellict  (^7^  (19U8)  gave  a solution  of  the  three- 
dimensional  field  equations  for  a transversely  isotropic  medium  in  terms 
of  two  stress  fbnctions,  each  satisfying  a second-orter  partial  differential 
eqoiation.  This  solution,  in  the  special  instance  of  rotational  stress 
syiatetiy  about  the  axis  of  elastic  symmetiy,  reduces  to  that  oidginated 
by  Lekhnitsky  Elliot's  approach,  which  is  evidently  incomplete  in 

the  absence  of  axisymmetzy  of  the  stress  field,  has  been  applied  by  him 
Q7],  and  by  Shield  Q9]]  to  the  solution  of  auiditional  significant 
special  problems  with  and  without  rotational  symmetiy.  Finally,  we  refer 
to  a recent  paper  by  Ana  Moisil  [[lO^  (1950),  in  which  a formal  scheme  due 
to  Gr.  C.  Moisil  [11]  is  used  to  extend  Galorkin's  |^63  stress-function 
approach  to  mediA  with  transverse  isotropy.  The  method  of  derivation 
employed  here  does  not  ass\ire  the  completeness  of  the  approach  which 

^See  [l],  p.  27U. 

^bi3  observation,  according  to  Mindlin  [|5],  is  due  to  H.  M.  Wester- 
gaard. 
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involves  three  stress  functions,  each  satisfying  a partial  differential 
equation  of  the  sixth  order.  No  attempt  is  made  to  effect  a reduction 
in  the  unnecessarily  high  order  of  the  equations  diaracterizing  the  stress 
functions,  which  is  an  essential  disadvantage  from  the  point  of  view  of 
applications,  and  the  case  of  axisymmetric  stress  fields  is  not  considered 
especially. 

It  is  the  purpose  of  the  present  paper  to  present  a systematic 
derivation  of,  and  to  supply  the  apparently  still  missing  conpleteness 
proof  for,  Lekhnitsky's  approach  to  problems  with  torsionless  axisynnetry 
in  a medixm  possessing  transverse  isotropy.  To  this  end  we  first  general- 
ize a theorem  of  Almansi  T12J  concerning  the  representation  of  the 
general  solution  of  the  equation  F = 0 in  terms  of  harmonic 

functions.  The  theorem  established  here  may  be  of  interest  beyond  the 
present  application.  The  material  in  Sections  3,  h,  although  primarily 
expository,  has  briefly  been  included  in  order  to  render  the  present 
paper  sensibly  self-contained. 
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2.  A Generalization  of  Almansl *s  Theorem 

In  this  section  we  prove  the  followinR  theorem: 


Let  R be  a region  of  the  (y3fZ)-plane  such  that  a stral^t  line 
arallel  to  the  z-axls  intersects  the  boundary  of  R in  at  most  two 


joints » Let  F^(p,z)  be  a solution  of 


in  X2 


(2.1) 


where 


= A+  cJdJ 

1 1 z 


A = 


(2.2) 


the  are  (not  necessarily  real)  constants,  and  the  functions  p.  are 


continuous  in  R with  Pj.  ^ 0,  Then  admits  the  representation^ 


_ 

Fn(^jz)  = * * I^(/>^z) 


(2.3) 


n - 1 


where t 


; J ^i^n-l  = °'  ^n^  =0' 

2 

and  m (0  ^ m ^ n-1)  is  the  number  of  the  coefficients  c^  (i 

2 

n-  1)  which  are  equal  to  c^. 

For  future  reference,  we  first  note  the  trivial  identities 

[a(z)B(;>,.)]  = Ar=B  ♦ c2(a^_^,B  ♦ 

, (i)  , , (i)  , (1) 

V “ '7^  =0. 


(2.h) 


= 1,  2, 


(2.5)- 


(2.6) 


Subscripts  preceded  by  a comma  denote  partial  differentiation  with 
respect  to  the  argument  indicated. 


li 


f* 

p 


and  observe  that  the  product  of  any  two  among  the  operators  D^,yV,  V^> 

2 

Vj  is  connutative. 

As  a further  preliminary,  we  show  that  if  P(jO,z)  satisfies  = 0 

and  m > 0 is  an  integer,  then 


vf(o'V) 


0 for  k = m - 1 

2®  [n  if  F for  k = a 


(2.7) 


Let  k = m - 1.  Clearly,  (2.7)  holds  for  m = 1.  We  proceed  by  induction. 
Thus,  assxne 


^2(m-l)  (jm-2  pj  ^ 


Then, 

= nn-2)z”"^  ■»  2z”"Vj=  0, 

by  hypothesis  and  (2.5).  Next,  let  k = ■•  Again,  (2.7)  evidently  holds 
for  m = 1.  Using  induction  once  more,  assume 

Then, 

= - l)z””^  F + 2z”“^  F^^J 

= If  F, 


by  hypothesis,  (2.5),  and  the  resxilt  for  k = m - 1.  This  corgsletes  the 
proof  of  (2.7). 
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Turning  to  the  proof  of  the  theorem,  we  need  to  show  that  there 
(n) 

exists  a function  F(y3,z)  which  satisfies 


= 0, 


n - 1 


1 r V?  = 0* 


(2.8) 

(2.9) 


llorooTer,  we  may  assume  without  loss  in  generality  that 

(i  = n - m,  n - m + 1,  ...,  n - 1), 
(i  " 1,  2,  ...,  n " m “ !)• 


2 2 
*=1  = ‘^n 


2/2 

®n 


J 


(2.10) 


In  view  of  (2.10)  we  may  write  (2.9)  as 
n - 1 


-m- 1 

m \ 

1 = 1 

7^  r.  ,(n)J 
^ 0 

(2.11) 


and  (2.11),  by  virtue  of  (2.6),  (2.7),  and  (2.8),  becomes, 

ttf  f'")  (c|  - c^).  (3.12) 

(n) 

Thus,  it  suffices  to  construct  a function  ¥{j>,z)  which  satisfies  (2.8) 
and  (2.12).  To  this  end  let 

- 1 


nr2n  ♦ m + 2 


T 


g(j3»z)  = 


i = 1 
n-m~  1 


Vi^n 


(2.13) 


cfTT  (4-cb 

* 11  a _ -1  1 n 


i = 1 


.-1 


the  operator  D being  defined  by 
z 


D"^G(jo,z)  = 


(2.11) 
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irtiere  is  the  value  of  z belonging  to  one  of  the  points  of  intex^ 
section  of  a straight  line  parallel  to  the  s-axis  with  the  boundary  of 
R*  Evidently,  = g satisfies  (2.12).  Furthermore, 

= 0 (2.15) 


because  of  (2.1).  Hence, 


2n-m-  3 


V^g  = 


(2a6) 


Now,  let  g^(j3fZ)  be  defined  by 


2n-m-3 

g,=  ^ 


'^2n-m-3  “ ^2n-m-3’ 


^2n-m-U  = ^2n-m-l,  ^<2.17) 


A\  = fjt  “ ®n  * 2)h^^2  = 0,l,...,2n-m-5), 


The  function  g^  so  determined  has  the  properties 


v;  *.  = v^g, 


(2.18) 


It  follows  from  (2.13),  (2.18)  that 


= e - 


g - g. 


(2.19) 


satisfies  both  (2.8)  and  (2.12).  This  completes  the  proof  of  the  theoreo. 

Successive  applications  of  the  foregoing  theorem  at  once  yield  a 

representation  of  any  solution  of 

n 

I r V?  F = 0 (2.20) 

i = 1 ^ 

in  terms  of  solutions  of  the  equations 


4 


8 


1 


2 .oVi)  _ Q 
2 


VI  T 


(i  — 1 j 2 j • • • » n) . 

In  particular,  if  all  Cy  (i  = 1,  2,  n)  are  distinct^ 

F = 


(2.21) 


n 


.(k) 


(2.22) 


•whereas  if  all  coalesce. 


n 


F = 


(2.23) 


In  connection  with  the  plane  problem  of  elasticity  theory  for  an 
orthotropic  medium,  special  interest  is  at-taelied  to  the  case  in  which 
(2.2)  assume  the  form. 


A=  ^ i 


^2  _ 3^  2 9^ 

i i?- 


(2.2U) 


Here  j3  and  z are  Cartesian  coordinates,  and  V^  is  a modified  two- 
dimensional  Laplacian  operator.  In  the  rotationally  symmetric  problem 
appropriate  to  a medium  with  transverse  isotropy,  on  the  other  hand,  we 
encounter  the  case  in  which 


A = — ? 

3y)  j) 


2 _ 1 a . ^2  3^ 

p 


(2.25) 


where  p and  z are  circular  cylindrical  coordinates  and 
modified  axi symmetric  Laplacian  operator. 


is  a 


Finally,  if  in  (2*20)  is  defined  as  in  (2.2U)  or  (2.25)>  and 


c^  = 1 (i  = 1,  2,  n),  then  (2.23)  becomes  identical  with  Almansi's 


^12]  representation  of  the  general  solution  of  = 0 in  terms  of 

,2 


hannonic  fimctions,  for  the  case  in  •which  V i®  •two-dimensional  or 
the  axisyrametric  three-dimensional  Laplacian  operator.  This  special  in- 
stance arises  in  the  well  known  applications  of  Almansi's  theorem  to  the 
plane  and  to  the  rotationally  sj’nunetric  problem  in  •the  classical  theory 
of  elasticity  of  isotropic  media. 
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3»  Stress-Strain  Relations  and  Strain-Energy  Fxmction  for  a Medium 
Tflth  Transverse  Isotropy 
Let  us,  temporarilj',  adopt  the  notation. 


^ ~ ^xx*  ^ ^ ^ ^ ~ ^ ~ 


± XX'  2 yy'  3 zz'  u 5 zx'  6 xy 


(3.1) 


Here,  T^t  •••  and  e^,  *.*  are  the  Cartesian  cooponents  of  stress 


and  "infinitesimal"  strain,  the  strains  being  defined  by 


®xx 


» • 


e _ = u + u , 

yz  y>2  z,y" 


(3.2) 


u , u , u 

L ^ y zj 


are  the  Cartesian  scalar  ccmponents  of  the  displace- 


•where 
ment  vector. 

The  general  linear  stress-strain  law  now  assumes  the  form. 


^ s e^, 


(3.3) 


in  which  i,J  range  over  the  integers  one  to  six,  and  the  usual  sxanmation 
convention  is  employed.  For  a homogeneous  medixnn  the  c.  . are  constants.^ 

A J 

A necessary  and  sufficient  condition  for  the  existence  of  a strain- 
energy  function  W(e^,  ...,  e^),  such  that 


. 'bH 
7i  - dej' 


(3.U) 


is  that 


= "ir 


Moreover, 


W = i c 


? ®i 


(3.5) 


(3.6) 


Note  that  the  present  definition  of  the  elastic  constants  c,  . 
coincides  with  Uial  adopted  by  Love  ^l]]^  p.  99. 

9 


TW 


10 


We  non  impose  the  condition  (transverse  isotropy)  that  the  stress-strain 
law  (3.3)*  (3*5)  be  invariant  under  an  arbitrary  rotation  about  the 
z-axis.  This  is  the  case  if  and  only  if 


c^j  = 0 (i  = 1,2,3;  j = (i  = U,5,C;  J i) 

°11  " °22»  °13  " °23*  " ^SS*  ^^66  " ^11  ~ ^12' 


K3.7) 


j 


With  the  notation. 


Cii  = °33  = ^6S  ~/^ 

c^2  = a - 


h (3.8) 


the  stress-strain  law  (3»3)»  (3*5) » subject  to  (3.7),  (3*8),  becomes 

+ (a  - 22?)  e„  + be„_ 

XX  XX  / yy  zz 

+ (a  - 2JI)  e^  be„„ 

' yy  yy  ^ ^ zz 


h (3.9) 


T-LZ  = V * 

^zr  ”/^®zr*  ^xy* 

whereas  (3*6)  appears  as 


2W  = a(ef^  + e^  ) + ae^  + 2(a  - 2JAj)e  e 

' XX  yy'  zz  xx  yy 

+ 2b(e + e )e  + /^(e^  + e^  ) -^ZT^  . (3.10) 

XX  yy'  zz  yx  zx.  xy 

The  stress-strain  relations  (3*9)  involve  five  independent  elastic  con- 
stants, of  which  ^ and  have  an  obvious  physical  meaning* 
Invertin(»  the  first  three  of  (3*9),  we  reach 
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Farthennore, 


E=5^ra  + 7^)-  (3.13) 

l^e  phsrsical  significance  of  the  four  new  elastic  constants  defined  in 
(3.12)  is  immediate  from  (3.11). 

A trivial  compiitation,  based  on  (3«10)»  yields  the  following  neces- 
sary and  sufficient  conditions  for  the  positive  definiteness  of  W: 

a > 0,  a > 0,  ^ > 0»  ^ >0,  aa  - - a^  > 0,  (3*lii) 

or,  equivalently, 

E > 0,  E > 0,^  > 0,^  >0,  -1<  >3<1,  i->)>2  (3.15) 

In  the  special  case  of  isotropy,  we  have 

a = 5 = 2/d  -»-A  = b = ;\  , 

=A  = ^(i  ♦ ';)')  • 


(3.16) 


With  reference  to  cylindrical  coordinates  (p,  0,  y),  where 


X = j3  cos  y,  y = yo  sin  y,  z = z,  (U.1) 

the  displacement  field  in  the  presence  of  torsionless  axisynmetry  about 
the  z-axis,  is  characterized  by 


a.2) 


^e  associated  field  of  strain  now  beccmes 


e = u , 

ff 


e = u -f  u . 
pZ  jO,Z  z,p* 


^yp  - ®zy 


and,  according  to  (3*9),  the  stress-strain  relations  are 


f'pf  = y ♦ '’V 

= a.yy  + (a  - Jj)e  + 

^ ■ (I1.I.) 

^zz  = * V 

The  stress  equations  of  eqviilibrium,  in  the  absence  of  body  forces,  by- 
virtue  of  (U.2),  {hm3)f  and  (U.U)  here  reduce  to 

/ - / 

/ 4 r'  — n=o 

pp,J>  yOZ,Z  O 

^ ^ ^ *a.5) 

/ 4 y"  4 ^ = 0 

OZ.O  'zz.z  o 


and  substitution  of  (li.2),  (h»3)>  into  (li.5)  yields 


ae 


*/^y,zz  + = o> 


- (/i  + b - i)u^_ 


zz 


= 0, 


in  nhich 


0 = e 


+ e. 


4 e = — 


1 d 


jf  n j,  dj> 


= : t 


_ , 4 U 

z,z 


KU.6) 


(li.7) 
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• Derivation  and  Ccwpleteness  of  Lekhnitsky*3  Stress-Function  Approach 
Siippose  u^{j)yz)  and  u^(y),z)  satisfy  the  equilibrium  equations  (U.6) 
in  a region  R of  the  meridional  half -plane  yj^O,  -oo<z  < oo,  and  let 
R have  the  property  that  straight  lines  parallel  to  the  coordinate  axes 
intersect  the  boundary  of  R in  at  most  ttro  points.  We  may  then  define 
two  ftinctions  U(y3,»)  and  V^,z)  by  aieans  of 


u = U , 


u = V 


(5.1) 


For  convenience,  we  introduce  the  operators 

^7*  = " ®z  = " I"  = “ ^ 

djO  papp  op  op 

D = — , 

z a 
dz 


(5.2) 


(5,3) 


being  the  axisymmetric  Laplacian.  Substitution  of  (5.1)  in  (li.6), 
with  the  aid  of  (1.7)  and  (5.2),  yields 

^ J^a  U + /XC^  + (/i  + b)  vj  = 0 

^ [(//  + b)  v|o  *yuvlt*  a>^]  = 0, 

or,  since  (5.1)  determine  U and  V only  within  arbitrary  additive 
functions  of  z and  p,  respectively, 

aV^  +//D^  + (//  + b)  D^  = 0 

^ + b)  V ^ + /UVlv  + = 0 

On  eliminating  first  V and  then  U between  (5»li),  we  reach 


(5.U) 


jtlu  = 0,  nv  = 0, 


2h 


(5.5) 
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the  operator  -TL  being  given  by 


* ^ " ® & 8 


(5.6) 


Consistent  with  the  notation  introduced  In  (2*2) > we  now  define  the 


operator  through 


v|  = (i  = 1,  2) 


(5.7) 


2 .2 


and  seek  to  determine  the  coefficients  c£,  in  such  a way  that 


rL=  v|. 


(5.8) 


To  this  end,  we  expand  the  right-hand  member  of  (5.8)  as  a polyncmial  in 
vl>  ^ equate  its  coefficients  to  the  corresponding  coefficients 
tqjpearing  in  (5*6) « This  leads  to 


_2  . .2  _ aa  - 2b/<.  - .2  .2  _ a 

c^  + C2  = ^ , c^Cg-  — , 

a/4 


(5.9) 


.2  2 


so  that  c£,  C2  are  the  roots  of  the  quadratic. 


Q(^)= 


4^  = 0. 


(5.10) 


The  inequalities  (3«lii)  imply  that  Q(^  ) cannot  have  a negative  or 
vanishing  root.  Conditions  (3«lU),  however,  do  not  preclude  the  exisb- 


2 2 

ence  of  complex  roots.  Indeed,  c£  and  Cg  are  real  and  distinct,  co- 


alescent,  or  conjugate  coioplex,  according  as 


aa  - (^  + 0* 


(5.n) 


In  the  special  instance  of  isotropy  we  conclude  from  (3«l6),  (5«10)  that 


c^  = Cg  = 1»  Tdience  (i  = 1»  2)  aixl  £1.=  V^. 


In  vienr  of  (5.5) j (5*8),  and  by  virtue  of  the  theorem  proved  in 


Section  2,  U admits  the  representation 

(1)  . (2) 

U = B + z B 


2 

B -=0, 


B 


(2) 


= 0 


(5.12) 


ufaere 


k = 


.0  if  ^ 

•1  if  c^=  c|. 


(5^3) 


(1)  (2) 

It  is  expedient  to  define  two  new  functions  k.y>,z),  k^,z)  by  means  of 


(1)  (1)  (2) 

B = k _ +kA.  , 


(1) 


V£  A =0, 


(2)  (2) 

® = \zz 

(2) 


A 0. 


(5.1U)‘ 


Fran  (5.1^),  (5.12),  (5.13),  (5.7)  follows 

p (1)  (2)-i 

U = [a  + z*^  a,,  J 


.2  _2  k ,,..2 


= + <=2  - 2>"2  * 


] 


(5.15) 


(5.16) 


Substituting  (5.15),  (5.16)  in  the  first  of  (5.U),  and  performing  two 
successive  z-integrations,  we  obt^lin 

r 2 2 k <2) 

^ + b)  V = I (ac£  -/^)A  + (ac^ -^)z*"  A^^ 


2 (2) 

- 2ac2kA  + z f . 


1 ♦ ^2]* 


(5.17) 


where  f^  and  are  functions  of  j>  alone.  We  first  treat  the  general 

case  in  which  the  degeneracy  ^ = - b is  ruled  out. 



The  existence  of  A'  and  A is  readily  verified  with  the  aid 
of  the  last  two  of  (5.12). 
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Case  1:  ^ + b 0 

In  order  to  determine  the  nature  of  the  functions  f^,  fg,  ire  substi- 
tute (5»16),  (5.17)  into  the  second  of  (5«ii).  This  leng^chy  computation, 
with  the  aid  of  (5.10),  (5.l3)j  (5«lU),  leads  to 


fi  = 0,  f j = 0 


Now  define  <J)(yj,z)  by 


. 1 r k / T 


From  (5.19),  (5.15),  (5.1.),  (5.1,  follows 


u = - (//  4-  b)  <b 

f ‘ 

= a (J)  + - a)  4^^ 


and  (5.13),  (5.U),  (5.18),  (5.19)  ijnply 


(5.18) 


J 


(5.20) 


rt<t)  =vlvl(p  =0. 


(5.21) 


We  have  shown  that,  barring  the  exceptional  case  yCC  + b = 0,  every 
solution  of  the  displacement  equations  of  equilibriun  (U.6)  admits  the 
representation  (5.20)  in  terms  of  the  stress  function  <J)  idiich  is  sub- 
ject to  (5.21).  Moreover,  a direct  substitution  confirms  that  every  dis- 
placement field  of  the  form  (5.20),  with  a particular  solution  of 

(5.21),  satisfies  (li.6)  even  if  yU.  + b = 0.  Equations  (5.20),  (5.21), 
except  for  an  m. essential  constant  multiplier  of  (p,  are  identical  with 
the  stress-function  approach  given  by  Lekhnitsky  ^3j,  and  independently 
arrived  at  by  Elliot  Applying  the  theorem  of  Section  2 to  (5.21), 
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and  recalling  (5«7)>  we  note  that  the  general  solution  of  (5»2l)  may  be 
written  as 


(f>(c,2)  = z/c^)  4- 


(5.22) 


where  (p^(p,z)f  (p^(j),z)  are  arbitrary  axisynmetric  harmonic  functions 
and  k is  given  by  (5.13) • The  cylindrical  components  of  stress  belonging 
to  the  displacements  (5«20)  are  obtained  with  the  aid  of  (li.3),  (U.U)  and 
may  be  omitted  here. 

In  particular,  if  the  medium  is  isotropic,  let 


1-2^ 


(5.23) 


In  view  of  (3.16),  Equations  (5.20),  (5.21)  here  reduce  to 


2/4 = 2(1  - X- 

= 0. 


Y (5.2U) 


which  is  Love's  stress-fxinction  approach  to  isotropic  problems  with 
torsionless  axisymmetry. 


Case  2;  yCl  + b = 0. 

In  this  special  instance  (5.U)  degenerate  into 
and  by  (5.10), 


(5.25) 


(5.26) 


fl 


See  H.  p.  276. 


Consequently,  according  to  (5»1)#  (5»7)»  the  general  solution  of  (U*6)  nor 


appears  as 


u = U , u = V 


P *P 

V^  = o, 


V^=  0 


M5.27) 


and  is  thus  again  e:q)ressible  in  terms  of  two  arbitrary  harmonic  functions. 
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